Spin Seebeck effect in antiferromagnets and compensated ferrimagnets 
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We theoretically investigate the spin Seebeck effect (SSE) in antiferromagnets and ferrimagnets, 
and show that the SSE vanishes in antiferromagnets but survives in ferrimagnets even at the magne- 
tization compensation point despite the absence of its saturation magnetization. The non-vanishing 
SSE in ferrimagnets stems from two non-degenerate magnons. We demonstrate that the magnitude 
of the SSE in ferrimagnets is unchanged across the magnetization compensation point. 
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I. INTRODUCTION 

Much attention is now focused on the thermal effects 
in spintronics, and the emergent research field of spin 
caloritronics is rapidly developing One of the most 
important issues in spin caloritronics is the spin Seebeck 
effect (SSE).— The SSE is the mechanism by which a 
spin voltage is generated from a temperature gradient 
in a magnetic material over a macroscopic scale of sev- 
eral millimeters^. Because the spin voltage is a potential 
for electron spins to drive spin currents, this spin voltage 
injects a pure spin current, i.e., a spin polarized current 
which is unaccompanied by a charge current, from the 
ferromagnet into an attached nonmagnetic metal. The 
inverse spin Hall effect (ISHE) 5,6 converts the injected 
spin current into a transverse electric voltage and hence 
the SSE is electrically detectable. Since its discovery in 
2008, this phenomenon has drawn much interest as a sim- 
ple way of generating pure spin currents that are needed 
for future spin-based technology;^ and the recent obser- 
vation of the giant SSE in InSb 9 has attracted a consid- 
erable attention. 

The SSE has been observed in various ferro- 
magnetic materials ranging from metallic ferromag- 
nets, NigiFeig^ and Co2MnSi, 10 to semiconduct- 
ing ferromagnet (Ga,Mn)As) 11 ' 12 to insulating mag- 
nets LaY 2 Fe 5 0i2^ and (Mn,Zn)Fe 2 4 ^ Although 
LaY2FesOi2 and (Mn,Zn)Fc204 are classified into fer- 
rimagnets in a rigorous terminology, the current under- 
standing of the SSE in these systems relies on a model- 
ing as ferromagnet a 1 5 i 16 because the low-energy magnetic 
properties relevant to the SSE are well described by a fer- 
romagnet modeling owing to the large gap between the 
acoustic and optical magnons. These observations have 
established the SSE as a universal aspect of ferromagnets. 

Besides ferromagnets, ferrimagnets and antiferromag- 
nets are known as prototypes of magnetic materials*^ A 
fcrrimagnct is an ordered spin system in which two sub- 
lattice magnetizations point in the opposite directions, 
and an antiferromagnet is classified as a special case of 
a ferrimagnet for which both sublattices have equal sat- 



uration magnetizations. Recently, there has been an on- 
going attempt to develop antiferromagnetic metal spin- 
tronics, and several experimental^ and theoretical^— 
work are already in progress. Regarding ferrimagnets, 
the intriguing characteristics of ferrimagnetic ordering 
are now drawing considerable attentio n 22 i 23 in develop- 
ing a ultrafast magnetization manipulation technique. 
Therefore, it is quite natural to ask whether the SSE 
can be observed in antiferromagnets and ferrimagnets. 

In this paper, we address the issue of observing the 
SSE in antiferromagnets and ferrimagnets. Especially, 
we focus on the SSE in ferrimagnets with magnetization 
compensation. A certain class of ferrimagnets are known 
to possess a magnetization compensation temperature 
Tm (angular-momentum compensation temperature Ta), 
at which the two sublattice magnetizations (spins) have 
the same magnitudes but opposite directions, leading to 
net-zero saturation magnetization (spin angular momen- 
tum)^— We show that two non-degenerate magnons 
give rise to the non-vanishing SSE at Tm or Ta despite 
the absence of net saturation magnetization or total spin. 
Also, we show that for a uniaxial antiferromagnet the 
SSE vanishes because the thermal spin injection by the 
two degenerate magnons is perfectly compensated. More- 
over, the SSE in an easy-plane antiferromagnet is shown 
to disappear because in this instance neither magnon car- 
ries spins. 

This paper is organized as follows. In Sec. II, we in- 
vestigate the SSE in uniaxial antiferromagnets as well as 
ferrimagnets with magnetization compensation. Next, in 
Sec. Ill we discuss the SSE in easy-plane antiferromag- 
nets. Finally, in Sec. IV we summarize and discuss our 
results. 



II. SPIN SEEBECK EFFECT IN UNIAXIAL 
ANTIFERROMAGNETS AND FERRIMAGNETS 

As a general model of ferrimagnets and antiferromag- 
nets, we consider the following Hamiltonian defined on a 
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FIG. 1: (Color online) Schematic view of a hybrid structure 
composed of a nonmagnetic metal (N) and a ferrimagnet (F) 
with two sublattices A and B. 

lattice composed of two sublattices A and B^ 

Hf = —J a ^ &aa ■ Saa' — Jb ^ Sbj ■ Sb.j' 

(i,i')eA (j,j')eB 

+ J AB Saa-Sbj + SHa + SHb, (1) 

(ieA,jeB) 

where J a and Jb (Jab) are the nearest-neighbor intra- 
sublattice (inter-sublattice) exchange integrals, and (, ) 
specifies nearest neighbor bonding (see Fig.QJ. The last 
two terms in Eq. (JTJ) for sublatticc L = A,B are given 
by SH L = Eiei^Mo^o ■ - • S Lti )% where 

fiQ is the Bohr magneton, Hq — —HqZ is the external 
magnetic field, g L and Dl are the effective g-factor and 
the anisotropy constant for sublattice L. 

First, we use the spin- wave approximation to diagonal- 
ize Eq. (J]). Following standard procedures^ using the 
linear Holstein-Primakoff transformation for spin opera- 
tors = S x L l ± iSl ti (L = A,B), the Hamiltonian © 
is diagonalized to be 

h f = {"t a \ a * + w « > ( 2 ) 
i 

where w± = \^{e A + ef ) 2 - 4^ ± ( e £ - c «), and the 

precise forms of e~, e q , and 77^ are given by e q = 
2z J A S A [l - 7g] + zqJabSb + (g A HoH + D A S A ) and 

= 2zoJ_B5B[l-7q] + 2; O^ABS'A + (-.g B A i O^O + -C , B'S'B)- 

Here, 7 q = Zq 1 ^2 s e lq s is defined by the sum over 
zq nearest neighbors of the original lattice, and rj q = 

JabVSaS b J2s' eiqS ' is defined by the sum over zq near- 
est neighbors of the sublattice A or B. In this paper, 
we assume a cubic lattice for simplicity. In Eq. ©, 
the magnon operators a q and (3 q are defined by the 
Bogoliubov transformation^ a q = u q a q + u~/3 q and 
bq = u q a q + u qPqi where and a q and 6 q are the 
Fourier transforms of operators a; = (2Sa)~ 5 i an d 
6 S ; = {2S B )-^S B l with 5a = |SU| and S B = \S B \, and 

,2 _2 ' 
U+ - U q = 1. 




FIG. 2: (Color online) Spin-wave spectra (Ho = 0) with q 
along the [111] direction calculated from Eq. ([2]) using pa- 
rameters for (a) a uniaxial antiferromagnet NiO, and (b) a 
compensated ferrimagnet Er3FesOi2. The wavevector q is 
measured in units of the inverse of the nearest-neighbor dis- 
tance. 

In Fig. the spin-wave spectra (Ho — 0) calculated 
from Eq. (|2|) for a uniaxial antiferromagnet NiO and a 
compensated ferrimagnet Er3FesOi2 are plotted. For 
NiO, we use Jab = 6.3 meV (J a = Jb = 0), D = 
0.1 meV, S A = S B = 0.92^- whereas for Er 3 Fe 5 0i 2 , 
we assign the net spin of the rare-earth ions (the fer- 
ric ions) to Sa (Sb) on a model cubic lattice, and we 
set J A = meV, J B = 0.68 meV, Jab = 0.19 meV, 
S A = 4.2, S B = 2.5, g A = 1.4, g B = 2.0 D A = 3.5 x 
10~ 3 meV, and Db — 3.0 x 10~ 4 meV to reproduce the 
Neel temperature TNeei = 556 K and the magnetization- 
compensation temperature Tm = 83 K i 28 ' 33 As is well 
known, the two antiferromagnetic magnons are degener- 
ate if Hq = 0, whereas the two ferrimagnetic magnons 
are non-degenerate because of the inequivalence of the 
two sublattices. 

We discuss now the SSE in uniaxial antiferromagncts 
and ferrimagnets modeled by Eq. ([1} . Note that a uniax- 
ial antiferromagnet can be modeled as a special case of 
a ferrimagnet. We consider a model shown in Fig. [TJ in 
which a ferrimagnet (F) and a nonmagnetic metal (N) 
are interacting weakly through the s-d exchange interac- 
tion at the interface. We assume that the ferrimagnet F 
has a local temperature Tp, and the nonmagnetic metal 
N has a local temperature Tjy. We analyze the SSE in 
the longitudinal configuration 3 ^ by employing the linear- 
response formulation of the SSE in a ferromagnet de- 
veloped in Refs. [IH and [35|. The s-d interaction at the 
interface is modeled by 

W-sd = ^2 (j^i (T i ■ Saa + J^&j ■ S B ,j^J(3) 

i,j£F/N-intcrf&cc 

where, for sublattice L = A.B. J^ d is the s-d exchange 
interaction at the F/N interface, er^ is the itinerant spin 
density operator in N. The total Hamiltonian of the 
system, H, is then given by 

H = H F +n N + H sd , (4) 
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FIG. 3: (Color online) Feynman diagram representing two 
processes relevant to the SSE in uniaxial antiferromagnets and 
ferrimagnets. (a) Spin current injected by a q magnons (I*). 
(b) Spin current injected by /3 q magnons The signs of 

/J" and 17 are opposite. The solid and wavy lines represent 
magnon and itinerant spin-density propagators, respectively. 

where Hat is the single- particle Hamiltonian of the con- 
duction electrons in N (see, e.g., Eq. (31) in Ref. I36I). 

The central quantity that characterizes the SSE is the 
spin current / s injected into TV, because it is proportional 
to the experimentally detectable electric field .Eishe via 
ISHE££ 




FIG. 4: (Color online) Temperature dependence of the SSE 
signal I s [red, Eq. Q], saturation magnetization M 3 [blue, 
Eq. l|9j)]. and total angular momentum Stot [green, Eq. (|10p ], 
calculated for a compensated ferrimagnet Er3FesOi2 using 
the same parameters as in Fig.[2](b). The case for a M s pinned 
by the anisotropy field is shown; the data is normalized by its 
value at T/Tw& e \ = 0.1. Inset: The case for a M s pinned by 
the external magnetic field is shown. 



-Eishe — dsnpJs x cr, 



(5) 



where #sh and p are respectively the spin-Hall angle and 
the resistivity of N, and J s = (I s /Ai nt )x is the spin- 
current density across the F/N interface having a con- 
tact area A int . Following Refs. [JJ| and [H, we calculate 
I s as the rate of change of the spin accumulation in N, 
i.e., I s — J2ieN(dt&i) where (• • • ) denotes the statisti- 
cal average. What is special in the present calculation 
is that we need to express the s-d interaction [Eq. 
in terms of the a q and /3 q operators ["±" branches in 
Eq. (J2J)], because these are the magnon operators in F. 
Following procedures presented in Appendix [AJ the spin 
current injected in N is expressed as 



2y/2 



^N N N F h , q 
+JJ(k,q)B?Ju> 



Re E / [^dCMM&H 



(6) 



where A K (B K ) is the Keldysh component of the inter- 
face correlation function between magnon operator a q 
(ffl) and the itinerant spin-density operator (see Ap- 
pendix [A]) , and we have introduced the shorthand no- 
tation J = ji. Here J^{k,q) is the effective s-d 
interaction written in terms of magnon operators, and 
the precise definition is given in the Appendix [Al 

We perform the perturbative approach in term of the 
s-d interaction at the interface to evaluate Eq. ((6|). Then, 
the spin current I s is given by the two diagrams shown 
in Fig. [31 and accordingly, I s has two terms: 



where if, representing the contribution from the ± 
branch, is expressed by 

it = ±E 8Arint H ( t 9)|21 / we*,*) 



k.q 



N N N F h 2 



xlwggiq,*) coth(^)-coth(^; 



.(8) 



(7) 



Here N lnt is the number of localized spins at the F/N 
interface, Nn (N f ) is the number of lattice sites in N 
(sublattice sites in F), {\J s f{k,q)\ 2 } = S A (J^f + 

Sb(J>^) 2 - In Eq- ©> X R {k,u) = Xn/(1 + X%k 2 - 
icjr s f) where ^jy, Aat, r s f are respectively the param- 
agnetic susceptibility, the spin-diffusion length, and the 
spin-flip relaxation time in TV, and Q±(q,ui) — l/(w — 
uj q + ia±uj) where a± is the damping parameter in F. 
Note that the signs of the spin current injected by the 
a q magnons (If) and that by the /3 q magnons are 
opposite. 

We first consider the SSE in a uniaxial antiferromag- 
net. As is depicted in Fig. [2] (a), the two magnons in a 
uniaxial antiferromagnet are degenerate if Hq = 0. More- 
over, owing to the equivalence of sublattices A and B. the 
s-d exchange interactions at the interface for these two 
sublattices are the same (J s ^ = Jjj). From these condi- 
tions we obtain |/+| = |/~| resulting in a null SSE due 
to Eq. ([7]), i.e., I s = 0. Thus, the SSE vanishes in a 
uniaxial antiferromagnet under a negligibly small exter- 
nal magnetic field because of the perfect compensation 
of the spin injection by the two degenerate magnons. 

We next consider the SSE in a ferrimagnet close to 
the magnetization compensation point, in which the two 
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magnons are no longer degenerate. Figure [4] shows the 
temperature dependence of the SSE signal I S (T) calcu- 
lated from Eqs. and flEJ for a compensated ferrimag- 
net Er3Fe50i2 by using the same parameters as in Fig. [5] 
(b). In Fig. [4] we also plot the saturation magnetization 

ieA " jeB 

to determine the magnetization compensation point de- 
fined by M s (Tm) = 0. In addition, we plot the total 
angular momentum 

Stat = (S Z ) (10) 

to determine the angular-momentum compensation point 
defined by St o t(T0 = 0. Here, S z is the z-component of 
the total spin S = Sa + Sb, i.e., 

ieA r jeB 

Clearly we see that the SSE signal is unchanged across 
both compensation points, either Tm ~ 0.15Tnc C i or 
Ta ~ 0.32Tn6c1- We performed the same calculation 
for several different choices of parameters, and confirmed 
that the SSE is unchanged across Tm and Ta. 



where e± = y/ (A% + 2B± ) {A± - 2B± ) , A± = 

2z JS cos 2d + g^ Q H sin 6 + \ D \ S T zqJS{cos 26 - 1)7,, 
B± = tz JS(cos26 - 1)7, - \D\S/2, and lq = 
Zq 1 ^2 s e lq S is defined by the sum over zq nearest neigh- 
bors. In the above equation, 8 is the canted angle of 
the sublattice magnetization, and the magnon operators 
and are defined by the Bogoliubov transformation 

^ (a, + 6_q) = u q i q + v q ^_ q and -±(a q - b_ q ) = 

x q(q + VqC-q) where u q - v q = 1 and x q - y q = 1 are 
real coefficients, and a q and b q are the Fourier transforms 
of Holstein-Primakoff operators and bi. As is seen in 
Fig. 3 of Ref. H3, the two magnons in the easy-plane an- 
tiferromagnet are not degenerate even when Hq = 0. 

Now we discuss the SSE in an easy-plane antiferro- 
magnet modeled by Eq. (|12[) by using the same proce- 
dure as in the previous section. As before, we consider 
a system in which an easy-plane antiferromagnet (eAF) 
having local temperature T e A_F and a nonmagnetic metal 
(N) having local temperature Tjv are interacting weakly 
through the s-d exchange interaction at the interface. In 
the absence of an external magnetic field, a direct calcu- 
lation shows that the spin current I s injected into N is 
identically zero, i.e., 

I. = 0. (14) 



III. SPIN SEEBECK EFFECT IN EASY-PLANE 
ANTIFERROMAGNETS 

In this section, we show that the SSE in easy-plane an- 
tiferromagnets vanishes under a zero magnetic field be- 
cause neither of magnons carries spins in easy-plane an- 
tiferromagnets. We consider the following Hamiltonian 
for easy-plane antiferromagnets^i 

HeAF = J S ?- S 3 

{ieA,jeB) 

D 



£ Y,^oH -S^--(z-Sf) 2 }(12 



L=A,B ieh 

where J is the nearest-neighbor exchange integrals, Hq = 
HqX is the external magnetic field, and g is the g-factor 
and D < is the anisotropy constant which selects the x- 
y plane as an easy plane. Note that the external magnetic 
field is applied along the x axis, and we assume Sa || z 
and Sb \\ —z when Hq = 0. Following Ref. we in- 
troduce the linear Holstein-Primakoff transformation by 
choosing the direction of each canted sublattice spin in 
the ground state as a spin quantizing axis. Performing a 
7r/4 rotation to the operators to separate the mixing of 
the two spin operators and using the Bogoliubov trans- 
formation, Eq. (fT2")) is diagonalized to be 



This is understood by investigating the z-component of 
the total spin S = Sa + Sb [Eq. (TITl) ]. In the case 
of a uniaxial antiferromagnet discussed in the previous 
section, the expectation value of S z is given by 



<s*> = (s A -jL £K«,}) - (s B - ^ £<fltfl 



N F 



(15) 

where a q and j3 q are the magnon operators defined in 
Sec. im Equation (|15l) means that the a q magnons car- 
ries spin one while /3, magnon carries spin minus one in 
a uniaxial antiferromagnet. On the other hand, the ex- 
pectation value of S z in an easy-plane antiferromagnet 
under discussion is calculated to be identically zero, i.e., 



(S z ) = 0. 



(16) 



Equation (1161) means that magnons in an easy-plane an- 
tiferromagnet are similar to a linearly-polarized photon 
and hence neither of magnons carries spins if Hq = 0. 

In the presence of a finite external magnetic field 
Hq = HqX with a nonzero canted angle 9, however, the 
x-component of S becomes nonzero. In this situation, 
the spin current I s injected into N is shown to have the 
polarization along the x axis, and its magnitude is given 
by 



eAF 



= h 



q 



a i<;i</ "F £ q CoC 



q SqS,q 



q Sq^Q 



(13) 



gyp H Q 
Jz 



+ 17) 



(17) 
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where 



N N N F h 2 



i. - J ui 



coth( 



k.q 

2/cbTjv 



) - coth( 



2kuT e AF ' 



(18) 
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Here, F±(q,u>) — l/(u> — e q + ia±ui) is the retarded com- 
ponent of the magnon propagator with a± is the damp- 
ing parameter. Note that the signs of 7+ and I~ are the 
same, in contrast to the case of a uniaxial antifcrromag- 
net. 

From Eqs. ([14]) and JT7]), we conclude that the SSE 
vanishes in an easy-plane antifcrromagnet if Hq = 0. 



IV. DISCUSSION AND CONCLUSION 

The main result of this paper is that the SSE in an- 
tifcrromagnets vanishes, whereas the SSE in ferrimag- 
nets persists and is insensitive to either magnetization or 
angular-momentum compensation effects. The interpre- 
tation is as follows. For the SSE to occur, the existence of 
the transverse fluctuations of the total spin, i.e., 5f ot and 
S% ot , is needed. For a ferrimagnet at Tm or Ta, fluctua- 
tions of S x ' y do not vanish even when Stot = or M s = 0. 
Therefore, ferrimagnetic magnons can always generate 
the SSE. Only for a uniaxial antiferromagnet, where the 
two magnons are degenerate, the SSE from the two de- 
generate magnons with the opposite sense compensates 
perfectly. Note that neither magnon in an easy-plane 
antiferromagnet carries spins. 

Our conclusion is not modified by considering the 
phonon-drag contribution to the SSE2& because, as dis- 
cussed in Refs. |H and 0, phonon drag can be taken 
into account by replacing Tp and Tn in Eq. ([HJ with an 
effective magnon temperature T F and effective spin accu- 
mulation temperature T^. We also note that the magnon 
excitations are well defined even at Ta- The presumed 
divergence of the magnon damping parameter at Ta— 
does not exist when we recall the condition justifying the 
expansion used in Rcf. 42;, where such an effect mani- 
fests itself as an enhancement of the damping parameter 
without any divergence (see Appendix [Bj . Note that the 
magnitude of magnon damping has less effect on the lon- 
gitudinal SSE, although it has a large influence on the 
transverse SSE. 

To conclude, we have theoretically investigated the 
SSE in antiferromagnets and ferrimagnets, and shown 
that the SSE vanishes in antiferromagnets whereas it 
persists at either the magnetization or the angular- 
momentum compensation points of ferrimagnets, despite 
the absence of its saturation magnetization or total spin. 
Because a fringing field by saturation magnetization is 
suppressed at the magnetization compensation point, 
this phenomenon can be useful for constructing a pure 
spin current device which is free from crosstalk of the 
fringing field. 



Appendix A: Linear-response expression of 
magnon-driven spin injection in ferrimagnets 

In this Appendix, we derive Eq. ([6]) in the main text. 
We consider a system described by the Hamiltonian (Q}, 
and calculate the spin current I s = J2 it£N (dt&i) ■ We 
use the momentum representation of <rf and calculate 
the quantity I s = \J N^s (d t (J^ o ) fc -s-o • The Heisenberg 
equation of motion for cr| gives 



d t cr 



fen 



where 



i£N/F(L) 




h.(cA,l) 



sublattice 



£ 

and we have used the relation [a%.,a , ] 
Taking the statistical average of the above quantity, the 
spin current is calculated to be 



= A,B, 



2 ± 
fiSfc k+k' 



- 4 ^ ReV 

y/N N N F h 



[j a i(k,q)A£ q (t,t) 



-J- A {k,q)B<(t,t) 



(A2) 



where J^(k,q) = J^(k,q)</SXu± 



J sd 



(k,q) 



Here, A< q (t,f) = -i(a q (t')a^(t)) and B< q (t,f 




— i(/3q(i') <T fe (*)) measure the interface correlation func- 
tions between the magnon operators (a q and f3 q ) and 
spin-density operator err . In the steady state the inter- 
face correlations A^ (t,t') and B^ q (t,t') depends only 
on the time difference t — t' . Introducing the frequency 



representation A< q (t,t') = <t A t q ^)^ (t ~ 
well as using the relationship A < = h [A K — A R 
we finally obtain Eq. ([5]) in the main text. 



-*') 
A* 



as 



Appendix B: Magnon damping near compensation 
points 

In this Appendix, we calculate temperature depen- 
dence of magnon damping parameter close to the com- 
pensation points and show that the magnon excitation 
is well defined even at compensation points without any 
divergences. We begin with two Landau Lifshitz Gilbert 
equations for sublattice L — A, 

dML - 1L M L x H, 



dt 



«*> Ml x^,(B1) 



M. 



s.L 



dt 
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FIG. 5: (Color online) Temperature dependence of the 
effective magnon damping parameter a e g [red, Eqs. (|B8|) 
and (|B9|) ]. saturation magnetization M s [blue, Eq. JSJ], and 
total angular momentum Stot [green, Eq. ()10[l ] . calculated for 
a compensated ferrimagnet Gd23Fe74.6Co3.4. The data is nor- 
malized by its value at T/Twee\ =0.1. 



where Ml is the sublattice magnetization with its mag- 
nitude given by Ml, Hl is the effective magnetic field, 
7l = gL^o/h is the gyromagnetic ratio, and oll is the 
Gilbert damping parameter. The effective fields are given 
by H A = Hq+H^-XMb and H B = H +H% n -XM A , 
where Ho = HqZ is external magnetic field, ff| n = H'^z 
and H B n = —H B n z are the anisotropy fields, and XMa 
and XMb are the inter-sublattice exchange fields with 
A = zoJab/{9a9b^'o)- Because we here focus on the 
uniform mode, the intra-sublattice exchange couplings 
X A = z J A /(g A ^o) 2 and X B = z Jb/(9b^o) 2 are dis- 
carded in Eq. (|BT|) . 

Below the magnetization compensation point we set 
M A = M A z + m,A and M B = —AI B z + ms, such that 
the effective fields can be written as 



H A = (H + Hf + XM B )z-Xm B , (B2) 
H B = (H - H% a - XM A )z - Xm A . (B3) 



fntroducing the representation E C A — —(Ho 
XM B ) and Ef = -(H - H%f - XM A ), and linearizing 
with respect to m A and m B , the Landau-Lifshitz-Gilbert 
equations are transformed to be 



dm A 
dt 



dm i 

dt 



= z x 



-y A (XM A m B - E c A m A ) + a A 



dm A 



= — z x 



j B (XM B m A + E c B m B ) + a B 



dt 

(B4) 
dm B ~ 



dt 



We introduce m ± = m x ± im y and substitute m \ (t) 



,+ „ — iujt 



into Eqs. (|B4| and (|B4| . Then we obtain 



(—iu> - a A uj + \~f A E c A )m A ~ - iX"f A MArrig = 0, 



(B6) 



(-iw + a B u + ij B E B s )m B ~ + iXj B M B m\ = 0. 

(B7) 

The eigenfrequency uj is determined by the equation: 

(lj - ia A u - j a E a s )(uj + ia B uj - -y B E e B ) 
+X 2 lAlB M A M B = 0. (B8) 



Above the magnetization compensation temperature, 
we set M A = —M A z + m A and M B = M B z + m B 
because we consider a situation in which the satura- 
tion magnetization is pinned by an external magnetic 
field. This situation can be analyzed by rewriting E A — 
-(Ho - HX - XM B ) and Ef = -(H + fff 1 + XM A ) as 
well as reversing the signs of a A and ag. We numerically 
solve Eq. (|B8|) by setting 



(B9) 



(B5) 



Figure [S] shows the temperature dependence of the 
effective Gilbert damping parameter a e g the lower fre- 
quency mode, calculated for a compensated ferrimagnet 
Gd 2 2Fe7oCo 8 .— We assign A sublattice for Gd ions 
and B sublattice for Fe ions, and neglect Co ions for 
simplicity. We set S A = 3.85, Sb — 3.5, g A — 2.0, 
g B = 2.05, H = 0.04 T, H^ = 0.0 T, H% n = 0.02 T, 
a A — 0.004, aB — 0.0039. The saturation magnetization 
and the total spin are calculated by using the mean field 
approximation by using J A b — 0.28 meV, J A — meV, 
and Jb = 0.34 meV to reproduece TnooI = 500 K. These 
parameters give Tm ~ 0.079Tn6c1 and T A w 0.09ITncc1- 
From the figure, we see that Gilbert damping constant is 
enhanced around the angular momentum compensation 
point Ta, but does not show any divergences. A small 
discontinuity at Tm stems from the fact that the spin 
quantization axis is reversed at Tm because of the pin- 
ning by the external magnetic field. The result obtained 
here justifies the statement in Sec. IIVI that the presumed 
divergence of the magnon damping parameter at Ta— 
does not exist when we recall the condition justifying the 
expansion used in Ref. 42, where such an effect mani- 
fests itself as an enhancement of the damping parameter 
without any divergence. 
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